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Abstract 

Dynamic stabilization by simultaneous primary and superharmonic resonances for high order nonlin¬ 
earity cancellation is demonstrated with an electrostatically-actuated, piezoresistively-transduced nanome¬ 
chanical resonator. We prove experimentally how the combination of both the third-order nonlinearity 
cancellation and simultaneous resonances can be used to linearly drive a nanocantilever up to very large 
amplitudes compared to fundamental limits like pull-in occurrence, opening the way towards resonators 
with high frequency stability for high-performance sensing or time reference. 
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Nanoelectromechanical systems (NEMS) are being developed for many applications such as 
ultrasensitive mas and gas sensing, subattonewton force detection at millikelvin temperature!^ 
or single-spin detectioiP. Recently and due to the early occurrence of their nonlinearitieP, reso¬ 
nant nano-scale mechanical devices have been used as platforms to explore fundamental questions 
in nonlinear classical mechanicP. 

Large amplitude vibrations are essential to the frequency stability of a resonator, and hence 
to the performance of a resonant nanosensor. The relative frequency noise spectral density of a 
resonator within its bandwidth is given b)^: 
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where Q is the quality factor, Sx is the displacement spectral density and Pq the displacement 
carrier power, i.e. the root mean square drive amplitude of the resonator 

The maximum amplitude is usually set by the onset of nonlinearit)!^, either for stability 
reason^ or noise mixing in the carrier side band^. Some work has been carried out in 
the past few years to overcome this limitation through hysteresis suppression by nonlinearity 
cancellatioiP^Ull jg usually achieved by adding a negative contribution to the third order 

(Duffing) nonlinear term by means of an electrostatic gate for instance. Unfortunately, this tech¬ 
nique is limited by higher order nonlinear terms which may be of importance even for small 
amplitudes, i.e. a fraction of the gap between the resonator and the electrode (typically below 20% 
of the actuation gap for electrostatically-actuated NEMSj^. The amplitude response can become 
a multivalued function of the frequency and five possible amplitudes for one given frequency is a 
clear signature of the physical significance of the fifth-order nonlinearitie^. It has been shown 
how the use of simultaneous primary (drive force at cOn, resonance at ujn) and superharmonic (drive 
force at a;„/2, resonance at Un due to the second order nonlinear term) resonances can tune and 
dynamically delay the onset of this multivalued, highly unstable behavioi^. 

In this paper, we demonstrate experimentally how the combination of these two techniques, 
namely the nonlinearity cancellation and simultaneous resonances, can be used to stabilize and 
linearly drive a nanomechanical resonator up to very large amplitudes compared to fundamental 
limits like pull-in occurrence. 

The NEMS device (see Figs, [^a) and[TJb)) consists of a cantilever beam of length I = 5 fxm, 
electrostatically driven and connected to two suspended piezoresistive gauges at a distance 
d = 0.15/ from its anchored end. The cantilever vibration induces stress in the piezoresistive 
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gauges which in turn is transduced into a resistance variation. It is fabricated on 200 mm Silicon- 
On-Insulator wafers using complementary metal oxide semiconductor-compatible materials and 
processes. The readout scheme and process is fully detailed irPI. 



FIG. 1: (a): Design parameters of the NEMS device (/ = 5 jim, a = 3.5 fim, g — 200 nm, d — 700 nm, 
h = 300 nm, Ig = 500 nm, hg = 80 nm). The gap g is the in plane distance between the cantilever 
and the electrode, (b): Food-coloured SEM image of the in-plane piezoresistive structure, (c): Electrical 
measurement scheme of the piezoresistive NEMS device for both 2fi down mixing technique (V^c = 0) and 
simultaneous resonances (Vdc ^ 0). Vac and Vdc are respectively the alternating current {AC) and the direct 
current {DC) voltages applied to the drive electrode. 14 is the bias voltage applied to the nano-gauges, at a 
frequency slightly shifted from that of 14c to down-mix the output voltage at a low frequency (typically a 
few lOkHz). PS, LPF are power splitter and low pass filter, respectively. 

The considered NEMS has shown excellent performance in terms of frequency stability, with 
more than IQQdB dynamic range, and achievable mass resolution of few tens of zeptograms 
when operating in the linear domairPHSl describe here how the electrostatic actuation of this 
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nanoscale device is suitable to implement the simultaneous resonances and increase its dynamic 
range even further. 

A variational approach, based on the extended Hamilton principleP^ has been used in order to 
derive the nonlinear equations of motion describing the flexural vibration of the device sketched in 
Figure[T];a). A reduced-order model is generated by modal decomposition transforming the contin¬ 
uous problem into a multi-degree-of-freedom system consisting in ordinary differential equations 
in timeP^. Assuming that the first mode is the dominant mode and neglecting the nonlinear damp¬ 
ing, the resonator dynamics can be described by the following dimensionless nonlinear ordinary 
differential equation for a single degree-of-freedom ai. 

= 4Kcl4cC cos (nt) + cos (2nt) (2) 

The dot denotes derivation with respect to time, cOn the natural frequency, 02 , 03 , a^, the coef¬ 
ficients of the quadratic, cubic and quintic nonlinearities respective!))^, fi the constant damping 
coefficient. The amplitudes of the first and the second harmonic of the drive frequency are re¬ 
spectively proportional to VacVdc and where 14c and 14c are the AC and DC voltages applied 
to the nanoresonator and C is a constant that depends on the geometric parameters. 

When il is tuned around the resonant response at a;„ is obtained by both the first and the 
second harmonic 212. More precisely, the 212-excitation mainly generates the primary resonance 
at iOn, while due to nonlinearities, the 12-excitation actuates a superharmonic resonance at a;„. The 
response at cOn is thus made of simultaneous primary and superharmonic resonances. The latter 
one is generated via a ’’slow” excitation compared to the resonant frequency. Experimentally, this 
can be achieved using the ”212 mode” readout scheme described in Fig. [TJc). By contrast, if 
the frequency doubler is removed and the drive frequency is tuned around cOn with 14 c 7 ^ 0 and 
Vac -C 14 c, an ”112 mode” is simply obtained for conventional primary resonance only. 

The device under test was wire-bonded to a radio frequency (RF) circuit board and loaded in 
a vacuum chamber for measurements at room temperature. In order to avoid signal shortage by 
parasitic impedances, a down-mixing techniquepl was used to readout the resistance variation of 
the two nanowire gauges at a lower frequency Ate. Differential read-out was performed with in- 
phase and out-of-phase bias voltages 14 (one gauge is under compressive stress while the second 
one is under tensile stress In order to analyze the dynamic behavior of the resonator under 
primary resonance and discuss the onset of nonlinearities, the conventional 112 mode (frequency 
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doubler removed) was first used. The cantilever was actuated using high DC voltages and a fixed 


AC bias voltage at the gauges (Vt, = 1.56 V peak-to-peak) so as to reach the nonlinear regime. 


The frequency response was measured using a lock-in amplifier in frequency sweep-up and sweep- 
down in order to obtain a full characterization of the resonator bifurcation topology. A constant 
quality factor Q = 5000 was measured on linear frequency resonance curves. 

Figure shows two nonlinear resonance peaks in IQ mode. The right-hand-side resonance 
curve was obtained for Vac = 150 mV and Vdc = 5 V. It displays a hysteretic softening behavior 
characterized by an amplitude jump-up at the bifurcation point B 2 and an amplitude jump-down 
at the bifurcation point B^, where the cantilever oscillation amplitude is around 75% of the gap. 

The left-hand-side resonance curve was obtained for Vac = 75 mV and Vdc = 8V. Since the 
actuation force is proportional to Vdc-Vac, both voltages have a direct consequence on the displace¬ 
ment amplitude: the resonance amplitude is lower than the one in the right-hand-side curve. The 
increase in the DC voltage also affects the resonator stiffness i.e. shifts down the resonance fre¬ 
quency and amplifies the nonlinear electrostatic stiffness as well. Under those drive conditions, 
the mechanical Duffing nonlinearity is negligible with respect to the electrostatic one resulting in 
an amplified softening behavior (a stiffening behavior was also observed under different exper¬ 
imental conditions). Despite a lower amplitude, the influence of the high-order nonlinear terms 
and specifically of the quintic term (05 oc V^,) is clearly visible: in frequency sweep-up, only one 
jump-up has been identified at the bifurcation point B 2 . However, in frequency sweep-down, two 
jumps have been observed: a jump-up at the bifurcation point Bi and a jump-down at the highest 
bifurcation point in the softening domain S 3 where the cantilever oscillation amplitude is around 
52 % of the gap. 

Obtaining a linear frequency response with this device at high amplitudes can be performed in 
principle under primary resonance only (10 mode) by cancelling out the cubic term {a^iVdc) = 0 ) 
in Equation Q and thus suppressing the hysteresis. 0:3 depends on the mechanical and electrostatic 
parameters and in particular on the DC drive voltage. Assuming that the second harmonic in 
Equation (|^ is negligible, the optimal DC drive voltage cancelling 0:3 can be computed fronP: 



(3) 


(The mechanical displacement between the anchors and the gauges being negligible, the device is 
dynamically equivalent to a resonant nanocantilever of length I — d) 
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FIG. 2: Nonlinear resonance frequency responses measured using a If] down-mixing technique and show¬ 
ing the location of the different bifurcation points {Bi, B 2 and S 3 }. Wmax is the cantilever displacement 
at its free end normalized by the gap. The hysteresis suppression is limited by high order nonlinear terms 
which leads to a highly unstable behavior for Vdc = 81/ and Vac = 75 mV. 

Applying Equation Q to the NEMS in use yields a relatively low optimal DC voltage 
^dcop 1 V. Reaching large amplitudes (larger than the standard nonlinearity offset) in a lin¬ 
ear fashion would then require significantly high AC voltages (14c > 0.5 V). Consequently, the 
assumption of negligible second harmonic terms is not valid anymore and Equation Q is not appli¬ 
cable. More importantly, the non-linear terms of order five influence the dynamic behavior of the 
nanomechanical resonator already at low amplitudes as illustrated by the left-hand-side resonance 
curve in Figure The nonlinearity cancellation is thus limited by potential dynamic instability 
such as pull-in phenomena. Finally, no linear and stable behavior could be observed using the 
cancellation of the third-order terms with the device under primary resonance only. 

In order to overcome this issue, we combined the use of this third-order term cancellation 
technique with the use of simultaneous primary and superharmonic resonances, which has been 
shown to delay the onset of the quintic non-linear term^. In this case, the amplitude X and phase 
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/? modulations of the response can be written as: 


X = - 

/3 = - 

”|”2(T — 
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where e is the small nondimensional bookkeeping parameter and a is the detuning parameter. 

In Equations Q and Q, both slow (superharmonic) and fast (primary) dynamics are present 
due to terms proportional to or « 2 , the coefficients of the quintic and quadratic nonlinearities. 
The AC voltage sets the 2Vt excitation while the DC voltage only amplifies the 0, excitation. 
Unlike the case of primary resonance only shown in Figure tuning the bifurcation topology in 
the present simultaneous resonances configuration can be obtained by simply increasing the DC 
voltage, i.e. amplifying the superharmonic excitation amplitude. 

In practice, simultaneous primary and superharmonic resonances can be implemented ex¬ 
tremely easily with the 2 r 2 downmixing readout scheme described in Fig [TJc), with Vdc 7 ^ 0. 
Moreover, high AC voltages should be used to increase the first harmonic cos{flt) in Equation 
Q. The contribution of the AC voltage in the nonlinear electrostatic stiffness was taken into 
account in Equation Q and the optimal drive DC voltage was analytically and numerically com¬ 
puted ias{Vdc,Vac) = 0) for Vac = 2U, resulting in Vdcop ~ 0.5 U. Fig. shows a resonance 
peak obtained with these values. In this case, the resonator behavior is almost perfectly linear up 
to extremely high amplitudes compared to the actuation gap (above 95% of the gap) and only a 
slightly softening behavior remains. Computations predict that for this set of parameters the non¬ 
linear electrostatic and mechanical stiffnesses are balanced and that the oscillation amplitude of 
the cantilever is close to 200 nm at its free end {i.e. the width of the gap). The observed remaining 
softening part under these computed drive conditions are attributed to fabrication uncertainties and 
approximations in Vdcop calculations. 

In this experiment, the second-order nonlinear term «2 giving rise to superharmonic resonance 
is used to control the stability of the nanocantilever around its primary resonance by retarding the 
effect of the fifth-order nonlinearity. The latter yields a highly unstable behavior with five possible 
amplitudes for one given frequency at amplitudes below 50% of the gap g as shown in Fig. 

We demonstrate here how the use of simultaneous primary and superharmonic resonances permits 
the dynamic stabilization of the nonlinear nanomechanical resonator in order to maintain a linear 
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FIG. 3: (a): Slightly softening frequency response measured with the 2f2 down-mixing technique at the 
optimal DC voltage V^c — 0.5 V for Vac — (I 4 = 1.56 V). (b): Zoom in the frequency region of the 

resonance peak. Unlike the peaks obtained under primary excitation only in Figure the combined use 
of third-order term cancellation with simultaneous primary and superharmonic resonances enable a stable 
and linear behavior of the nanomechanical resonator at very large amplitudes close to the actuation gap (the 
onset of unstability being here above 95 % of the gap g). 

behavior at large-amplitude vibrations close to the gap. At this vibration level, the free end of the 
cantilever reaches a distance of a few nm to the electrode without a damageable pull-in as would 
occur without the use of simultaneous resonance. 

To experimentally confirm such an oscillation level, the DC voltage was gradually increased 
from 0.5 U up to 2 V. Obviously, the vibration amplitude of the cantilever should increase, eventu¬ 
ally hitting the actuation electrode. Besides, the cantilever behavior should subsequently become 
more softening. Figure shows the frequency response obtained for V^c = 2 V: the increase in 
the electrostatic softening nonlinear stiffness is evidenced by the frequency shift between the two 
bifurcation points (softening domain) which is significantly enlarged compared to Fig. How¬ 
ever the output signal at the peak is still Kut — 5.4mU, i.e. the same value as at Vdc = 0.5 V. 
The observed plateau is not due to any kind of electrical saturation and is very reproducible across 
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FIG. 4: Frequency response with Vdc — Vac — 2V (14 = 1.561^). The horizontal branch intercepting the 
bifurcation point B 3 is not due to any saturation in the readout scheme. It rather shows that the cantilever 
free end comes in physical contact with the actuation electrode in a reproducible way. This further confirms 
the high vibration amplitude reached by the device in Fig. 

multiple measurements. We attribute this to the fact that the resonator reached a stable vibration 
amplitude equal to the gap (200nm). At one given bias voltage, this can be interpreted as a mea¬ 
surement of the experimental output voltage sensitivity per unit of displacement of the cantilever 
free end and confirms the very high amplitude reached in the measurement shown Fig. [^compared 
to the gap. 

In this letter, the nonlinearity cancellation has been demonstrated on a nanomechanical res¬ 
onator electrostatically actuated based on piezoresistive detection (160 nm thick) fabricated using 
a hybrid e-beam/deep ultraviolet lithography technique. The device has been characterized using 
a in down-mixing technique in its nonlinear regime to investigate the stability domain of the can¬ 
tilever dynamics under primary resonance. It has been shown that the high-order nonlinearities 
limit drastically the operating domain of the nonlinearity cancellation. In order to overcome this 
limit, the sensor has been actuated under its primary and superharmonic resonances simultaneously 
by using a readout scheme. Specifically, the resonator can be linearly operated for displace- 
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ments well beyond the multistability limit and almost up to the gap, while retarding undesirable 
behaviors and suppressing pull-in phenomenon. 

This technique opens promising perspectives for time-reference or sensing purposes: the fre¬ 
quency stability of the device should be significantly improved compared to its operation in linear 
regime below the bistability limit. In particular, the mass resolution one can expect with this de¬ 
vice without reaching any damaging amplitude and this technique is about 20 zeptograms (five 
times smaller than the performance reported iiP). Specific devices will be fabricated in the near 
future to substantiate this fact. 
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